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Age-dating a star-burst with GEMINI/CIRPASS observations of 
the core of M83 
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Abstract. We present preliminary results from a set of near-IR integral field spectroscopic observations of the cen- 
tral, star-burst, regions of the barred spiral galaxy M83, obtained with CIRPASS on Gemini-S. We present maps in the 
Pa/3 and [Fell] 1.257 /im emission lines which appear surprisingly different. We outline the procedure in which we will 
use Pa/3 emission line strengths and measures of CO absorption to determine the relative and absolute ages of individual 
star-forming knots in the central kpc region of M83. 
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1. Introduction 



The central star-burst in M83 (NGC 5236) has been stud- 
ied photometrically by Harris etal. (2001) using HSTAVFPC 
images in the broad-band near-UV and optical, as well as 
narrow-band Ha and H/3to derive colours and line equiva- 
lent widths for 45 clusters. Despite the excellent spatial res- 
olution of these observations, optical photometric analyses 
such as these suffer from patchy dust extinction, a redden- 
ing vector which parallels the evolutionary tracks in the two- 
colour diagram, and selection effects which tend to exclude 
the very youngest (f <5 Myr) clusters which have strong emis- 
sion lines, but only a weak stellar continuum. Additionally, it 
is not possible, on the basis of broad-band colours alone, to 
distinguish an instantaneous burst of star formation from a 
constant star formation rate. 

To help overcome these effects, we use the Integral Field 
Unit (IFU) of the Cambridge Infra-Red Panoramic Survey 
Spectrograph (CIRPASS; Parry etal. 2000)' to measure both 
the Pa/3 equivalent width, as well as the CO(6,3) spectro- 
scopic index over the entire star-forming arc and nucleus of 
M83. As demonstrated by Ryder etal. (2001) for MlOO, the 
combination of two such diagnostics allows one to constrain 
not just the age, but also the burst duration, for each cluster 
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2. The CIRPASS spectrograph 

CIRPASS is a fibre-fed IR J+H band spectrograph with 
Multi-Object (MOS) and Integral Field (IFU) modes. The 
cryogenic camera is liquid nitrogen cooled while the spec- 
trograph optics are housed in a refrigerated cold room cooled 
to -42° C. The fibre-fed design introduces considerable flexi- 
bility of operation for CIRPASS, allowing observations on a 
range of telescopes (from the 8 inch laboratory test telescope 
used for system verification, the 4m class Anglo-Australian 
and William Hershel telescopes to the 8m Gemini south tele- 
scope) with minimal modification to the focal plane interface. 

In MOS mode 150 independent fibres of ^L6arcsec 
diameter are available over a field-of-view (FoV) of 
17 X 17arcmin to 40 x 40arcmin (fibre aperture and FoV are 
dependent on the details of the telescope focal station occu- 
pied). Fibres are typically deployed in pairs to target up to 75 
objects at a time while simultaneously sampling the sky. The 
telescope is then nodded between the two fibres of each pair 
to allow sky subtraction. 

CIRPASS performs integral field spectroscopy using an 
array of 490 macrolenses which make up the Integral Field 
Unit (IFU), shown in Figure |2l A macrolens array is used 
in favour of a microlens array to reduces losses from Fo- 
cal Ratio Degradation (FRD, Parry etal. 2000). Two broadly 
rectangular fields of view are available, 13x4.7arcsec and 
9.3 X 3.5 arcsec. The respective lens scales are 0.36 arcsec and 
0.25 arcsec. Finer scales are possible if taking advantage of a 
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Fig. 2. Here we see the Integral Field Unit (IFU) of CIRPASS. 
This unit is attached at the focal plane of the telescope. The 
broadly rectangular array of fibres can be seen through the 
field lens. 

During H band observations, marginally increased variabil- 
ity was observed in OH sky emission features, resulting 
in sub-optimal sky subtraction. This is, however, not be- 
lieved to compromise the observational dataset. Three over- 
lapping regions, close to the optical nucleus of M83, were 
observed (see Figure The 0.36 arcsec lens scale was 
used giving a FoV per pointing of ^^13x4. 7 arcsec. Obser- 
vations where obtained in both the J and H bands centred 
at Accn — 12530 A and Accn — 15590 A. The wavelength 
ranges cover the Pa/3 and [Fell] 1 .257 emission lines in J 
and the [Fell] 1.644/im line and CO(6,3) molecular absorp- 
tion feature in the H band. Two 900 sec background limited 
exposures were taken in each band at each location, inter- 
leaved with 900 sec off-target observations for sky subtrac- 
tion. Observations are performed using a Non-Destructive 
Read (NDR) mode. Three sequences of 10 detector reads 
were done at 5 min intervals throughout the integration, to 
facilitate cosmic ray rejection and to suppress read-noise. A 
small dither of ^0.36 arcsec (one IFU lens) was performed 
between on-target frames to allow the removal of the small 
number of detector defects. 

3.1. Data Reduction 

The data have been processed using the CIRPASS IRAF pack- 
age^. A data reduction cookbook detailing the package tasks 
is under development^. Data reduction proceeds in the usual 
manner. A reset frame correction, analogous to bias subtrac- 
tion for a CCD frame, is performed. Cosmic ray rejection is 
carried out for the 900 sec frames using the IRAF IMCOM- 
BINE task and CCDCLIP pixel rejection to combine the time 
series data obtained during each NDR observation. A pixel- 
to-pixel flat field correction and a bad pixel mask are derived 



Fig. 1. HST colour composite image from Harris etal. 2001 
(F300W, F547M and F814W) overlaid with the outline of 
the CIRPASS IFU indicating the three telescope pointings 
targeted. The optical galaxy nucleus is the bright spheroidal 
knot to the North-East of the centre of the image, at the end of 
the central IFU pointing. A naive interpretation of this three 
colour image suggests that much of the star formation activity 
has taken place in the southern arc of the nucleus, with iso- 
lated pockets of activity elsewhere. Using IFU spectroscopy 
observations one can survey the entire nuclear region, looking 
for the signatures of star-formation, past and present, without 
the bias inherent in previous long slit studies which have, by 
necessity, only targeted regions highlighted by broad band 
imaging techniques. 

telescope equipped with an adaptive optics system. In IFU 
mode CIRPASS is typically operated in medium resolution 
mode (/?^5000) giving a dispersion of ^2.2 A/pixel. Such a 
resolution allows CIRPASS to capitalise on the intrinsically 
dark sky background between the OH air-glow lines (Martini 
& DePoy 2000), which are responsible for much of the high 
sky brightness observed in the near-IR 

Prior to September 2003, an Hawaii IK HgCdTe array 
has been used in the CIRPASS camera, giving a single shot 
spectral range of ^2200 A at a two pixel resolution of 4.4A. 
The detector was upgraded to an Hawaii 2K array in Septem- 
ber 2003 and has been used successfully for science during 
CIRPASS MOS mode observations at the Anglo-Austi-alian 
Telescope (AAT) during October 2003. The larger detector 
provides complete coverage of the J or H bands, between 
the atmospheric water absorption troughs, in single observa- 
tions. Alternatively, many programms may benefit from in- 
creased resolution without the loss of wavelength coverage 
that would be inherent with the IK array. 

3. Observations of the nuclear region of M83 

During March 8'^ and 9*'' 2003, observations were under- 
taken of M83 (NGC 5236), with CIRPASS in IFU mode 
at the cassegrain focus of the 8m Gemini-South telescope. 
Both nights were photometric with sub-arcsecond seeing. 



^ The CIRPASS IRAF package can be obtained from 
http://www.ast.cam.ac.uk/~optics/cirpass/docs.html 

^ A data reduction guide for the CIRPASS IRAF can be found at 
.http://www.ast.cam.ac.uk/~o ptics/cirpass/datared/cookbook_snl987a.php 
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from non-dispersed frames with the detector illuminated by 
a low power tungsten lamp. Interleaved sky frames are then 
subtracted from target frames, this step removing the need 
for the subtraction of a dark frame. The spectra are extracted 
using the optimal extraction routine described by Johnson 
efa/. (2003). Spectra are also extracted from a tungsten il- 
luminated dome flat field frame. The 2D spectral images are 
then divided by the extracted 2D flat field frame to account for 
fibre-to-fibre variation. A wavelength solution is derived from 
the the co-added sky frames using the OH air-glow lines us- 
ing the line list of Maihara etal. (1993). The spectra are then 
transformed to a common wavelength scale using the IRAF 
LONGSLIT package tasks FITCOORDS and TRANSFORM and 
pairs of dither positions are mosaiced. Flux calibration and 
removal of telluric absorption features is performed using the 
spectrum of an F2V star, HD 128299, observed at an airmass 
similar to the target observations. The stellar spectrum is fit- 
ted by the appropriate black body spectrum (rcff=6750K). 
Data analysis is performed within an interactive software 
suite written in the IDL programming environment. 



3.2. Image registration 

The / and H band data were recorded on different nights. In 
order to register the data a continuum image is created, recon- 
structing the IFU image across the full wavelength range in 
each band. Primary features in the final continuum image are 
the southern star-forming arc and the optical nucleus. Using 
these features we can register the data to within the individual 
IFU lens resolution (^^0.36 arcsec). 



3.3. Variance array propagation 

The two dimensional spectra are critically sampled on the de- 
tector with a FWHM of '--^ 1.8 pixels. The peak-to-peak sepa- 
ration of adjacent spectra is ?; = 2.0 pixels. In order to ex- 
tract spectra with such close separations, leading to a sig- 
nificant overlap between the wings of adjacent spectral pro- 
files, an optimal extraction technique is implemented solv- 
ing for the maximum likelihood light distribution of each fi- 
bre triplet (a fibre and its two nearest neighbours). The cen- 
tral position of each fibre, relative to a set of observations 
of well spaced {y > 100 pixels) calibration fibres observed 
at each instrument setting prior to the science exposures, is 
accurately known from extensive laboratory tests. 

This process (and that of accurate emission line profile 
fitting) requires the propagation of the expected variance as- 
sociated with each pixel of the detector With 900 sec ex- 
posures, the noise in CIRPASS observations is dominated 
by the count rate from thermal background emission below 
1 .67 fj,m. A short pass blocking filter prevents thermal emis- 
sion long- ward of 1.67 /xm entering the camera optics. At the 
small number of wavelengths affected by strong OH air-glow 
lines, noise is significantly enhanced by Poisson noise from 
the high intensity lines. 



3.4. Equivalent width measurements 

We fit profiles to the three prominent emission lines: Pa/3, 
[Fell] 1.257 /im and [Fell] 1.644/im using a reduced 
technique and assuming a single Gaussian profile. Such a 
simple profile is believed to be appropriate at the level of spa- 
tial and spectral resolution. During line fitting a spectral mask 
is used to weight the fit against pixels affected by strong OH 
air-glow line residuals. This OH avoidance technique is a fun- 
damental source of the sensitivity of CIRPASS. 

OrigHa, Moorwood & Oliva (1993) describe the CO(6,3) 
absorption band and provide a prescription for measuring the 
equivalent width of the feature. They define the local contin- 
uum using a linear fit to two continuum points which straddle 
the CO band. The equivalent width is then integrated over 
a predefined pass band. The width of this band (1.6175 /im- 
1.6220/im, 45 A or 20 pixels) is of the order of the velocity 
structure seen in the nuclear regions of M83 (±150kms~^, 
±8 A, ±3.6 pixels). It is therefore not desirable to simply 
define a single systemic velocity-adjusted pass band for the 
entire observational data set. To overcome this we define a 
smooth background velocity field based on the observed ve- 
locity of the Pa/3 emission line. 

4. Summary of results 

Quantitative analysis is still in progress, but already we see 
evidence for a much more complex star formation history 
than that hinted at by the single long-slit near-IR spectrum 
ofThatteeffl/. (2000). 

Figure |3l shows two example spectra. The left hand 
spectrum is taken from a region visible in narrow band 
HST/NICMOS Paa imaging with extended emission but 
which does not show signs of active star formation loci in 
broad band WFPC2/NICMOS imaging. This highlights one 
of the key goals of these observations, to perform an accu- 
rate census of the star formation in the nuclear region of 
M83, without bias towards bright continuum regions, such 
as those targeted by Harris efa/. (2001). In contrast to the 
strong Pa/3 emission, little CO absorption is seen in the H 
band spectrum of this region (not shown in Figure O, sug- 
gesting a young (< 5 Myr) star formation region. 

The right hand spectrum of Figure |3l is taken from the 
optical nucleus of M83. Strong CO absorption is visible in 
this H band spectrum which, when taken in comparison with 
a very weak Pa/3 emission line strength in the J band spec- 
trum (not shown), is suggestive of an old stellar population 
(> 15 Myr). By comparing the equivalent width ratio for 
Pa/3 and CO(6,3) we will determine the age of star-formation 
across the nucleus of M83. 

Figure |4] compares the distribution of Pa/3 line flux, with 
that of the [Fell] 1 .644 /xm. The former provides a measure 
of the massive star formation rate within the last ^5 Myr, 
while the latter is thought to be a tracer of shocks from 
supernova explosions within the past 10"* years (Alonso- 
Hererro etal. 2003). The [Fell] emission is rather clumpy, 
and confined mainly to the periphery of the more diffuse 
Pa/3 emission, suggesting that massive star formation is being 
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Fig. 3. Two example spectra are shown. 

Left : In the / band we see the strong Pa/3 emission, a signature of hot young stars and active star formation, along with the 
[Fell] 1.257 fim line. This spectrum is taken from a knot of strong hydrogen emission. While the knot is visible as diffuse 
emission in HST/NICMOS Paa imaging, it is not evident in broad band continuum observations at other wavelengths. 
Right : The CO(6,3) molecular absorption features are clearly visible in this H band spectrum taken from the bright optical 
nucleus of the galaxy (lower spectrum). The location of the centres of three of the CO band heads, after correction for redshift, 
are marked by dotted vertical lines. Two dashed vertical lines indicate the extent of the CO(6,3) band as defined by Origlia, 
Moorwood and Oliva (1993). A polynomial fit to the continuum level is also marked. The upper spectrum shows the OH 
air-glow emission spectrum from which a mask is created to weight spectral fitting algorithms. The CO features are believed 
to arise in the atmospheres of cool G/K and M red giant stars and can be used, in combination with the Pa/3 emission, to 
constrain the age and burst duration of star formation in this region. 
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Fig. 4. A comparison of the Pa/3 line flux (left), with that of the [Fell] 1.257 /im line (right) is made above. See discussion in 
Sectional 



propagated into regions of undisturbed gas by the passage of 
the supernova blasts themselves. 
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